Position-Space Description of the Cosmic Microwave Background 
and Its Temperature Correlation Function 
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We suggest that the cosmic microwave background (CMB) temperature correlation function C(9) as 
a function of angle provides a direct connection between experimental data and the fundamental cos- 
mological quantities. The evolution of inhomogeneities in the prerecombination universe is studied 
using their Green's functions in position space. We find that a primordial adiabatic point pertur- 
bation propagates as a sharp-edged spherical acoustic wave. Density singularities at its wavefronts 
create a feature in the CMB correlation function distinguished by a dip at 8 ~ 1.2°. Characteristics 
of the feature are sensitive to the values of cosmological parameters, in particular to the total and 
the baryon densities. 



The cosmic microwave background (CMB) radiation 
provides the best probe today of the early universe and 
a number of fundamental astrophysical constants [jlJ-Q ■ 
Dynamical evolution of primordial perturbations mani- 
fests itself in the form of "acoustic peaks" in the CMB 
temperature angular power spectrum C/. After over 
a decade of studies, the physical content of the peaks 
has become qualitatively understood |p|-pT|. Neverthe- 
less, one has to rely on standard numerical codes [ ^2|Jl3| 
to establish a quantitative connection between the CMB 
anisotropy pattern and cosmological parameters. This 
connection is particularly important now that high- 
precision CMB measurements have become reality. 

In this letter we show that previously unnoticed 
interesting phenomena are unraveled by considering 
the radiation-matter dynamics in position rather than 
Fourier space. Traditional CMB anisotropy analy- 
ses begin with the Fourier expansion of spatial inho- 
mogeneities and studying time evolution of individual 
Fourier modes Jl6|-^8[. The position space approach, 
while offering a formally equivalent description, differs in 
its physical interpretation, associated calculational meth- 
ods, and its implications for data analysis. As we show 
here, a new side of CMB physics is revealed in position 
space. One goal of our work is a deeper physical un- 
derstanding of both the CMB power spectrum C/ and 
the angular correlation function C(6). We demonstrate 
that acoustic evolution of perturbations before recombi- 
nation, as viewed in position space, produces sharp fea- 
tures in C (9) . These signatures not only provide a simple 
interpretation of the acoustic peaks, but they may also 
enable fast and accurate extraction of cosmological pa- 
rameters directly from C{9). 

As a simplified model for the essential physical pro- 
cesses, we consider the evolution of potential and den- 
sity perturbations in gravitationally interacting photon- 
baryon and cold dark matter fluids. We use the con- 
formal Newtonian gauge jl8],[l9) to describe gravity and 
assume adiabatic (isentropic) primordial fluctuations as 
preferred by the present data H|Lq] and the simplest in- 



flationary models. The photons are assumed to be tightly 
coupled to electrons and baryons by Thomson scatter- 
ing until recombination at redshift z rcc ~ 1100. Neu- 
trinos are treated like photons. The disregard of neu- 
trino free-streaming and photon diffusion introduces a 
I0%-20% error |o|,|lJ for CMB temperature anisotropy. 
This model is helpful for developing understanding. Our 
numerical results for the correlation function and the an- 
gular spectrum are obtained a full calculation with CMB- 
FAST m including all relevant effects. 

Long before recombination during the radiation era, 
the growing mode of gravitational potential pertur- 
bation is described in Fourier space by </>(k, r) = 
3 [sin(q)/<7 3 — cos(q)/<7 2 ] </>o(k), q = kc s T, where r is con- 
formal time, k is the wave vector corresponding to a 
comoving coordinate r, and c s — 1/Vo is the speed of 
sound in the radiation era |l6| . A delta function primor- 
dial perturbation <p(r,T — ► 0) = S^ 3 \r) is represented in 
/c-space by constant initial amplitude <fo(k) = 1. Fourier 
transforming the resulting 0(k, t) we obtain the three- 
dimensional Green's function for the gravitational poten- 
tial in the radiation dominated epoch, 



^ 3} (r,T) = ^(c s r)- 3 9(c s T-r), 

47T 



(1) 



where 6(x) = {I, x > 0; otherwise}. Eq. (Q) illus- 
trates the essential property of the position-space Green's 
functions — the perturbation is localized within the acous- 
tic sphere, with a sharp edge at the sound horizon 
r = S(r)^J T c s (r')dr'. 

When cold dark matter (CDM) and baryons are taken 
into account, we split the total gravitational potential 
4> = <p r + 4>c into a part connected to the photon-baryon 
plasma, <f> r> and a part due to the CDM, <fi c . The sources 
of the potentials <f> r and <f> c are given respectively by the 
radiation-baryon and CDM energy densities relative to 
the net zero-momentum hypersurfaces, the analogues of 
Bardeen's Jll| gauge-invariant density variable e m . The 
potentials <p r and <j) c unambiguously determine the val- 
ues of all the density and velocity perturbations in our 
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model j22j . Their evolution is described by a pair of cou- 
pled equations, 

4>r + E 4= r,c( A «^ + B ri<t>i) = <? s V^r , (2) 
4>c + T, l= r,c( Aci & + = 

where Aij and Bij are simple rational functions of r p2[ . 

The first of eqs. (Q) shows that a primordial isentropic 
perturbation at a point propagates as a spherical acous- 
tic wave in the photon-baryon plasma with the sound 
speed c 2 s = d Pl /d{ Pl + p h ) = (l/3)/[l + (3p b )/(4p 7 )], 
where Pl and pb are the photon and baryon mean energy 
densities. As the sound wave passes by, its gravitational 
effect perturbs the CDM causing it to evolve as described 
by the second of eqs. (^|). Since perturbations propagate 
differently in the radiation and the CDM components, 
inhomogeneities do not remain isentropic for r > . 

We integrate eqs. (^|) numerically along their char- 
acteristics starting from the radiation-era solution (|l|) 
at T; n i t ~ 10 -4 T rec , and increasing the space-time 
grid spacing with time in proportion to the acoustic 
horizon size S(t). It is convenient to solve eqs. (||) 
in only one spatial dimension with the initial condi- 
tions (jy^ix) —> S(x) and vanishing entropy pertur- 
bation as t — > 0. The resulting one-dimensional, 
or plane-parallel, Green's function for the total poten- 
tial, cjr- x > (x, t), is related to the conventional transfer 
function (j>{k,T) as 4>(k,r) = J <j>^ (x, r) e~ lhx dx. The 
potential Green's function in three dimensions equals 
0(3)(r ;T ) = -(2Trr)- 1 d(t> (1) (r,r) /dr. A superscript "(1)" 
will be used for any quantity (e.g. density perturbation 
or temperature anisotropy) obtained from the isentropic 
initial condition (f>(x,T — ► 0) = 5(x), and such quanti- 
ties will similarly be referred to as "Green's functions." 
The Fourier transform of any Green's function gives the 
corresponding fc-space transfer function. 

Fig. 1 shows the result for an effective temperature 
Green's function A^g (x) describing the combined intrin- 
sic and gravitational redshift contributions to CMB tem- 
perature anisotropy. Within our tight coupling assump- 
tion, it contains delta function singularities at the wave 
fronts given by (3/4) [1 + (3p b )/(4p 7 )]" 1/4 S (\x\ - S), 
where S (r) = c s (T')dT' . The contribution of the delta 

function to J A^ (x, r rec ) dx exceeds the absolute value 
of the contribution from the regular part by ~ 40%, 
showing that the singularities play a major role in CMB 
anisotropy. 

The massive baryons coupled to the radiation fluid re- 
duce the sound speed that decreases the acoustic horizon 
at recombination. They also drag radiation out of the 
x = region in Fig. 1 being repelled by a kink in CDM 
potential <pc at x = 0. This effect causes the prominent 
"gully" in the middle of the transfer function in Fig. 1. 

The peculiar features of the position space transfer 
functions are manifested in the angular correlation func- 
tion 
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FIG. 1. The position space transfer function for the in- 
trinsic (jS-y) plus gravitational redshift contributions to the 
CMB temperature anisotropy AT/T at recombination, for 
fi m = 0.35, Qa = 0.65, and h = 0.65. The vertical spikes rep- 
resent the wave-front delta function singularities described in 
the text. Baryon drag effects are evident from comparing the 
solid and dashed lines. 



C(9) ee (AT(n 1 )AT(n 2 )) = ]T 



(2Z + 1) 



4-7T 



GPi{cQs0), (3) 



where AT(nj) is the temperature anisotropy in the di- 
rection ni, and 9 is the angle between fii and n 2 . The 
connection between C (9) and the Green's functions may 
be illustrated in the large angle limit when AT(hi)/T w 
A ff(ni, T rcc ) = ~8 7 + 4>, 8 7 = Spj/pry, for 9 exceed- 
ing a few degrees. Writing each AT(fij) as the con- 
volution of the Green's function in Fig. 1 with a pri- 
mordial potential field cf>o, assuming that recombination 
occurs instantaneously, and averaging over the primor- 
dial fluctuations with power spectrum (</>o(ki)</>o(k 2 )) oc 
(5( 3 ) (ki +k 2 )fc" 3 - 4 , we get 



(4) 



where K is an analytically known function of x\ + a; 2 , 
angle 9, the comoving distance to CMB photosphere R, 
and the spectral index n s J22|. The exact formula for 
C(9) applicable to subdegree scales should also include 
Doppler and integrated Sachs- Wolfe effects and integra- 
tion along the line of sight Eq. neglects these 
contributions but is sufficient to give a good qualitative 
description of the effects considered below. 

From eq. (Q), the delta function singularities pro- 
duce a singular C(9) behavior in the vicinity of 9 S = 
2arcsin(5/i?), where S = 5(r roc ) ~ 140 Mpc is the 
acoustic horizon size at recombination. For a flat uni- 
verse and the scale- invariant power spectrum n s = 1, 
C{9)-C{9 S ) ~ -Asln|As|, As = sm(9/2)-S/R. When 
the two observed points come into acoustic contact at the 
critical angle 9 S , dC(9 s )/d9 = +oo in the instantaneous 
recombination approximation. At this angular separa- 
tion, a spherical sound wave emitted at r = from the 
point halfway between the observed points on the pho- 
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FIG. 2. Angular CMB temperature correlation function. 
The parameters of the ACDM model are f2 m = 0.35, 
S7a = 0.65, h = 0.65, n 3 — 1. For this set of parameters 
and the displayed Q^,h 2 values, 9 S ~ 1.17° and 1.22°. The 
thin solid curve is calculated using eq. (4). All the other plots 
are obtained in full CMBFAST calculation, except for the 
dotted line that artificially cuts the dip for its Fourier analy- 
sis in Fig. 4. The whole feature in C{9) is shifted to smaller 
angles in the open model (dash-dotted line). Baryon density 
variation modifies the slope for 9 < 9 s /2 (dashed vs. solid 
line) . 

tosphere has, by recombination, just reached the two ob- 
served points separated by a distance 2S. 

Paradoxically, the initial acoustic contact anti- 
correlates the temperature fluctuations causing a dip in 
the correlation function just below 9 S . As evident from 
Fig. 2 , the dip is also present in the full result obtained 
from a CMBFAST calculation (wide solid line) although 
smoothed by the Silk damping and the finite thick- 
ness of the recombination photosphere. An anomaly in 
the correlation function at this angular separation was 
pointed out earlier in [^4| and J2J|. At even smaller an- 
gles the correlation function rises steeply. The dip occurs 
when the comoving separation of the observed points on 
the photosphere is about 2S ~ 280 Mpc, providing a 
"standard ruler" in the early universe. The comoving dis- 
tance to the photosphere R and thus the observed angle 
9 S are determined by the rate of the subsequent Hubble 
expansion and by the cosmic curvature, which both de- 
pend on the total density parameter f2. The dash-dotted 
line in Fig. 2 demonstrates the shift in the visible location 
of the dip for an open model. 

What is the physical origin of the dip in the correla- 
tion function at the acoustic contact? According to our 
earlier discussion, A^(r) = — (27rr)~ 1 9A2j/<?r, which 
contains the derivative of a delta function. When con- 
volved with the primordial potential field </>q, S'(r — S) 
relates AT(r) to the gradient of 4>o(r). The vectorial na- 



ture of V</>o results in contributions of opposite signs to 
the temperature of two points located at distance S in 
opposite directions. 

A larger primordial power spectrum n s leads to a 
greater magnitude of all small-scale features in the cor- 
relation function. In the instantaneous recombination 
approximation, C(9) — C(6 S ) varies as As/(1 — n s ) for 
n s < 1 and As/|As|™ s_1 for n s > 1. Thus, the dip at 
9 = 9 S is enhanced for n s > 1 and suppressed for n s < 1. 

For 9 < 9 s /2 the correlation function slope becomes 
sensitive to the baryon density . This can be seen from 
comparing the dashed and the solid lines in Fig. 2. At 
such a small angular separation, the Doppler and ISW ef- 
fects provide significant contribution to the temperature 
correlation. Nevertheless, the fib dependence is largely 
determined by the same intrinsic plus gravitational red- 
shift contributions shown in Fig. 1. It is the interference 
of the baryon induced "gully" in the middle of one of the 
transfer functions in eq. (jij) and the delta function con- 
tributions of the other that causes the steeper C(9) rise 
at 9 < 9 s /2. The effect is almost linearly proportional 
to Cl^h 2 and may be useful for measuring the baryon 
density value using CMB data. 

The finite extent of the position-space Green's func- 
tions results in oscillations in their Fourier transforms, 
the fc-space transfer functions. These oscillations appear 
as the famous acoustic peaks in the CMB angular spec- 
trum Ci . The characteristic period of the oscillations 
is A; = ttR/S ~ 300 for f2 m = 0.35, n A = 0.65, and 
h = 0.65. The wavefront and x = Green's function 
singularities in Fig. 1 do not map directly to a particular 
distinct feature in C; spectrum because they are com- 
pletely delocalized by Fourier transformation. 

An interesting connection between features in C (9) and 
the acoustic peaks in C; is established by artificially cut- 
ting out the dip in C{9) and replacing it by a smooth 
interpolation with no local minimum as shown by the 
dotted line in Fig. 2. Fig. 3. shows that this operation 




FIG. 3. Angular power spectra for the ACDM models of 
Fig. 2. The dotted line gives Ci corresponding to the correla- 
tion function of the first model with the dip cut out. 
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eliminates all the acoustic peaks besides the first one. 
The first peak is closely connected with the sharp rise 
in C(9) at 9 < 9 S . Increasing the baryon density we 
increase the magnitude of the slope of C (9) . The magni- 
tude of the first peak grows accordingly (dashed vs. solid 
lines in Fig. 3) leading to the known reduction in the 
ratio of second to first acoustic peaks. 

In conclusion, the Green's function method provides a 
promising alternative to the conventional Fourier analysis 
of perturbations in the early universe. It offers new ana- 
lytical methods and intuition for the familiar phenomena 
such as the acoustic oscillations and the dependence of 
peak ratios on fl^h 2 . It also reveals new features that 
are localized and easily noticeable in position space but 
which were stretched over a multitude of wavenumbers in 
the Fourier domain. In particular, the singularities of the 
Green's functions lead to potentially measurable features 
in the CMB angular correlation function. 

Transfer functions are computed faster in position 
space than in Fourier space because the integration effort 
is concentrated within compact regions of space interior 
to the acoustic horizon. Computation of A^(x, T lec ) in 
Fig. 1 takes less than one second for 0.01% numerical ac- 
curacy on a PC. We are working to extend the method to 
integration of the Boltzmann equation for free-streaming 
neutrinos and for photons after recombination. 

Experimental data are collected in position space. As 
we have attempted to demonstrate, position space is also 
more appropriate for describing the acoustic dynamics 
of CMB radiation. It may thus seem natural to remain 
in position space when comparing observational results 
with theoretical predictions. The utility of the correla- 
tion function C(9) for data analysis has been recently 
considered by Szapudi et al. p3). Data analysis using 
the correlation function is complicated by the covariance 
of observational estimates of C{9) |27|]. Such correlations 
are eliminated for estimates of C\ made with uniformly- 
sampled all-sky maps. Nonetheless, the correlation func- 
tion may still be a valuable experimental tool because 
it concentrates much of the cosmologically essential in- 
formation into localized features. This may enable data 
processing efforts, and observing strategy, to be focused 
on the most relevant angular features. 
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